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Abstract. The spectrum of the self-adjoint Schrodinger operator associated 
with the Kronig-Penney model on the half-line has a band-gap structure: its 
absolutely continuous spectrum consists of intervals (bands) separated by gaps. 
We show that if one changes strengths of interactions or locations of interaction 
centers by adding an oscillating and slowly decaying sequence which resembles 
the classical Wigner-von Neumann potential, then this structure of the abso- 
lutely continuous spectrum is preserved. At the same time in each spectral 
band precisely two critical points appear. At these points "instable" embedded 
eigenvalues may exist. We obtain locations of the critical points and discuss 
for each of them the possibility of an embedded eigenvalue to appear. We also 
show that the spectrum in gaps remains pure point. 



1. Introduction 

In the classical paper |vNW29j von Neumann and Wigner studied the one- 
dimensional Schrodinger operator with the potential of the form '^'^'"(^^'^^) g^Yid dis- 
covered that such an operator may have an eigenvalue at the point of the continuous 
spectrum A = w^. Since then such potentials permanently attracted the interest of 
different authors [SMI lMi73l iRSTl iBiOTl iBiMl IHKS911 iKNOTl iNOTl iLTOl [JSTO] . 
Potentials of the type '^'^'"j^'^^) ^jth 7 £ (0,1) also possess such a property and 
they are often called Wigner-von Neumann potentials. 

Consider the operator on R+ that formally corresponds to the differential 
expression 

riGN 

with the boundary condition 

(1.1) ip{0)cos>c-ip'{0)sm>c^O, 

at the origin, where ao € M, d > 0, >f € [0, tt) and dx is the delta-distribution 
supported on the point x E R+ (for the mathematically rigorous definition of the 
operator see Section [^TTl cf. |AGHH05l Chapter III]). The operator Hj^ is self- 
adjoint in L^(R+). It describes the behavior of a free non-relativistic charged 
quantum particle interacting with the lattice dN. The constant ao characterizes 
the strength of interaction between the free charged quantum particle and each 
interaction center in the lattice. The spectrum of the operator has a band-gap 
structure: it consists of infinitely many bands of the purely absolutely continuous 
spectrum and outside these bands the spectrum of is discrete. The investigation 
of such operators was initiated in the classical paper jKP31] by Kronig and Penney. 

In the present paper we study what happens with the spectrum of the Kronig- 
Penney model in the case of perturbation of strengths or positions of interactions 
by a slowly decaying oscillating sequence resembling the Wigner-von Neumann 
potential. Let constants c?, ao, c, w, 7 and a real- valued sequence {<7n},ieN be such 
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that 

(1.2) c? > 0,Q!o G M,a; e ^0, — ^ ,7 e 1 ,{qn}net'S 

Model I: Wigner-von Neumann amplitude perturbation. We add a discrete Wigner- 
von Neumann potential to the constant sequence of interaction strengths. Let the 
discrete set X = {x„ : n e N} supporting the interaction centers be defined as 

(1.3) Xn ■= nd, n G N, 

and let the sequence of interaction strengths a — {an}neti be as below: 

/-, ~ csin(2a;n) 

(1.4) a„:=ao + ^ + Qn, neN. 

We study the self-adjoint operator ^ ~ which formally corresponds to the dif- 
ferential expression 

neN 

with the boundary condition (jl.ip . This operator reflects an amplitude perturbation 
of the Kronig-Penney model. 

Model II: Wigner-von Neumann positional perturbation. We change the distances 
between interaction centers in a "Wigner-von Neumann" way, i.e. we add a sequence 
of the form of Wigner-von Neumann potential to the coordinates of interaction 
centers leaving the strengths constant. Let the discrete set X — {S;„ : n € N} be as 
below: 

/-, ^\ ^ , csin(2a;n) 
1.5 Xn:=nd+ ^- - + qn, n € N, 

and the sequence of strengths be defined as 

(1.6) a„ :— ckQ, n gN. 

We study the operator ^ g which formally corresponds to the differential ex- 
pression 

neN 

with the boundary condition (jl.ip . This operator reflects a positional perturbation 
of the Kronig-Penney model and describes properties of one-dimensional crystals 
with global defects. We also mention that local defects in the Kronig-Penney are 
discussed in |AGHH05l §IIL2.6]; situations of random perturbations of positions 
were recently considered in [HITlOj . 

The essential spectrum of the operator H^^ has a band-gap structure similar to 
the case of Schrodinger operator with regular periodic potential: 

oo 

gUH^) = ^UH.) = U ([A+ -1, A^- _i] U [A2-„, A+ ]), 

where A+ < Af < A^ < A+ < • • • < A+ _i < \^^_^ < < ^tn < ■ ■ ■ 

The locations of boundary points of the spectral bands are determined by the 
parameters ao and d. Namely, the values A^ are the n-th roots of the corresponding 
Kronig-Penney equations 

Ls iVx] = ±1, 



where 

s\n{kd) 



(1.7) Ls{k) ~ cos{kd) + a^- 
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Ls{k) 



1 

cos(a;) 

— cos(w) 
-1 




Figure 1. The curve is the graph of Ls; bold intervals are bands 
of the absolutely continuous spectrum; crosses denote the critical 
points {A±„}„gN- 

For the details the reader is referred to the monograph [AGHHOSl Chapter III]. 

In the present paper we show that the absolutely continuous spectra of the 
operators ~ and ^ ~ coincide with the absolutely continuous spectra of the 

non-perturbed operator H^. However the spectrum in bands may not remain purely 
absolutely continuous. Namely, at certain (critical) points embedded eigenvalues 
may appear. In each band there are two such points. The critical points A^^j. in 
the n-th spectral band are the n-th roots of the equations 

= ± cosw. 

The illustration is given on the Figure I. 

For the operators considered we give exact conditions which ensure that a given 
critical point is indeed an embedded eigenvalue for some x g [0,7r). This can 
occur only for one value of >f G [0, tt). The possibility of the appearance of an 
embedded eigenvalue at certain critical point depends on the rate of the decay of the 
subordinate generalized eigenvector. We calculate the asymptotics of generalized 
eigenvectors for all values of the spectral parameter A including the critical points, 
except the endpoints of the bands. We also show that the spectrum in gaps remains 
pure point. 

Our results are close to the results for one-dimensional Schrodinger operator with 
the Wigner-von Neumann potential and a periodic background potential. Such 
operators were considered very recently in |KN07[ INSlli IKSll) . 

To study spectra in bands we make a discretization of the spectral equations and 
further we perform an asymptotic integration of obtained discrete linear system us- 
ing Benzaid-Lutz-type theorems ,BL87] . As the next step we apply a modification 
of Gilbert- Pearson subordinacy theory |GP87) . To study spectra in gaps we ap- 
ply compact perturbation argument in Weyl function method, see, e.g., |DM95] . 
|B(4Pn8) . 

The reader can trace some analogies of our case with Jacobi matrices. The 
coefficient matrix of the discrete linear system that appears in our analysis has a 
form similar to the transfer-matrix for some Jacobi matrix and the Weyl function 
that appears in our analysis takes Jacobi matrices as its values. 



4 



VLADIMIR LOTOREICHIK AND SERGEY SIMONOV 



The body of the paper contains two parts: the prehminary part — which consists 
of mostly known material — and the main part, where we obtain new results. In 
the preliminary part we give a rigorous definition of one-dimensional Schrodinger 
operators with (^-interactions (Section I2.ip . show how to reduce the spectral equa- 
tions for these operators to discrete linear systems in (Section 12. 2p . provide a 
formulation of the subordinacy theory analogue for the operators considered (Sec- 
tion Further we formulate few results from asymptotic integration theory for 
discrete linear systems (Section I2.4p and provide a characterization of the essential 
spectra of considered operators using the Weyl function associated with certain 
ordinary boundary triple (Section I2.5p . In the main part, in Section l3.ll we study 
a special class of discrete linear systems in and find asymptotics of solutions 
of these systems. After certain technical preliminary calculations in Section 13.21 
we proceed to Section [?751 where we obtain asymptotics of generalized eigenvectors 
for Schrodinger operators with point interactions subject to Model I and Model II. 
Further we pass to the conclusions about the spectra in bands putting an emphasis 
on critical points. Section r3.4l is devoted to the analysis of spectra in gaps by means 
of the Weyl function. 

Notations. By small letters with integer subindices, e.g. we denote sequences 
of complex numbers. By small letters with integer subindices and arrows above, 
e.g. itn, we denote sequences of C^-vectors. By capital letters with integer indices, 
e.g. Rn, we denote sequences of 2 x 2 matrices with complex entries. We use 
notations £p(N), ^p(N, C^) and £p(N, C^^^) foj. spaces of summablc (p = 1), square- 
summable (p = 2) and bounded {p = oo) sequences of complex values, complex 
two-dimensional vectors and complex 2x2 matrices, respectively. We denote for 
a self-adjoint operator its pure point, absolutely continuous, singular continuous, 
essential and discrete spectra by Upp, fiac, fsc, foss and fid, respectively. We write 
T e 6oo in the case T is a compact operator. Other notations are assumed to be 
clear for the reader. 

2. Preliminaries 

2.1. Definition of operators with point interactions. In this section we give 
a rigorous definition of operators with i5-interactions, see, e.g., |GK85[ lKo89] . Let 
a = {a„}ngN be a sequence of real numbers and let X = {xn : n £ N} be a discrete 
set on R-|_ ordered as < xi < X2 < ■ ■ ■ ■ Assume that the set X satisfies 

(2.1) inf \xn+i — Xn\ > and sup \xn+i — Xn\ < oo. 

nGN „gN 

Denote also xq := 0. In order to define the operator corresponding to the expression 

riGN 

and the boundary condition 

-0(0) cosx — ?/''(0) sinx — 0, 
consider the following set of functions: 

^-^^ ■■= --^^^'^ ^^>oc(K+ \ x\ , . e n| 

and let the operator H3^^x,a be defined by its action 

H^,x,c.4' := -4'" 

on the domain 

domiJ^,x,a := {V- e L'^{R+) n Sx,a- ip" G L^{M.+ ),^p{0)cos}< = ^P'{0)smx} . 
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The spectral equation Tx^ai^ = A?/; is understood as the equation —tp"{x) — X'ip{x) 
for E Sx,a- The latter is equivalent to the following system: 

- ^j"{x) ^X^{x), a; G K+ \ X, 

(2-2) 

1p{Xn + ) = IpiXn-), ■lp'{Xn + ) ij/ {Xn-) + an1p{Xn-), n € N. 

The equation (|2.2p has two linear independent solutions which are called generalized 
eigenvectors, li ip G satisfies (|2.2p and the boundary condition, then ip is 

an eigenfunction of H^^x,a- 

2.2. Reduction of the eigenfunction equation to a discrete linear system. 

In this subsection we recall rather well-known way of reduction of the spectral 
equation (|2.2p to a discrete linear system. Let the discrete set X and the sequence 
of strengths a be as in the previous section. To make our formulas more compact 
we introduce the following notations 

Sn{k) := sin(fc(a;„+i - Xn)) and c„(fc) := cos(fc(x„+i - a;„)), n e Nq. 

For a solution -0 of (|2.2p corresponding to A = fc^ we introduce the sequence below 

^„:=ip{xn), neNo. 

If the condition (EH]) is satisfied, then by jAGHHOSi Chapter III, Theorem 2.1.5] 
for k E C+ U R and n E N such that s„_i(fc), s„(fc) 7^ it holds that 

(2-3) -k( ^ + ^] + L + + ^) V„ = 0. 

ySn-l[k) S.a(k) J y Vs„_i(fc) Sn(k)J J 

Inversely, solutions of the eigenfunction equation on each of the intervals (a;„, Xn+i) 
can be recovered from their values at the endpoints x„ and Xn+i'- 

, . _ sin(fc(x„+i - x)) + sin(fc(2: - X„)) ^ ^ r . ^ n 

^ ' ^^^>~ sin(fc(x„+i ~ x„)) ' ^eFn,a:„+iJ. 

The reader may confer with |E97| . where a more general case of a quantum graph 
is considered. 

Instead of working with recurrence relation (|2.3p we will consider a discrete linear 
system in C^. Define 

it — ( ^"-1' 

V Sn 

Then (|2.3D is equivalent to 

(2.5) it n+l ^ Tn{k)lt n 

with 

/ 1 \ 

(2.6) Tn{k) := s„(fc) sin(fc(a;„+i-:c„_i)) a„s„(fc) 

\ s„-i(fc) s„-i(fe) ^ k J 

The coefficient matrix of this system Tn{k) is called the transfer- matrix. 

2.3. Subordinacy. The subordinacy theory as suggested in |GP87| by D. Gilbert 
and D. Pearson produced a strong influence on the spectral theory of one-dimensional 
Schrodinger operators. Later on the subordinacy theory was translated to differ- 
ence equations [KP92j . For Schrodinger operators with i5-interactions there exists a 
modification of the subordinacy theory, see, e.g., (SCS94] which relates the spectral 
properties of the operator H^^x.a with the asymptotic behavior of the solutions of 
the spectral equation (|2.2p . Analogously, to the classical definition of the subordi- 
nacy |GP87| we SRy tliRt Si, solution ipi of the equcitioii Tx^a 

ip — A?A is subordinate 
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if and only if for any other solution ip2 of the same equation not proportional to ipi 
the following limit property holds: 

We will use the following propositions to find location of the absolutely continuous 
spectrum. 

Proposition 2.1. |SCS94[ Proposition 7] Let H^^x.a be the self-adjoint operator 
corresponding to the discrete set X and the sequence of strengths a as in Section \2.1\ 
Assume that for all A C (a, h) there is no subordinate solution for the spectral 
equation Tx,a'4' = Then [a, &] C (j{Hy^^x,a) o,nd u{H^^x,a) is purely absolutely 
continuous in {a,b). 

Proposition 2.2. Let the discrete set X = {a;„ : n € N} on R+ he ordered as Q < 
xi < X2 < ■ ■ ■ ■ Assume that it satisfies conditions Ax, := inf„gN |x„+i — Xn\ > 
and Ax* :— sup„gf^ |x„+i — x„| < oo. Let a ~ {ctn}nefi be a sequence of real 
numbers and let A G M. Assume that liminf„_j.oo |sri('\/A)| > 0. // every solution 
of the equation Tx.aip — ^see (|2.2p ) is bounded, then for such A there exists no 
subordinate solution. 

Proof. Let ipi be an arbitrary solution of Tx.ai' — AV'- Differentiating (|2.4p one 
gets in the case s„(fc) ^ 

K(^)l< '^'('^"'+Jf'+^'^ forxe[x„,x„+i]. 

There exists N{k) such that inf„>^(;j) |sn(fc)| > 0. Since the sequence {^n}neN 
is bounded, one has that tp'i{x) is also bounded for x > Xx{k)- Obviously for 
X ^ XN(k) it is bounded too, since it is piecewise continuous with finite jumps at 
the points of discontinuity. Let ip2 be any other solution of Tx.a''P = Ai/i, which is 
linear independent with -01 . It follows that there exists a constant C > such that 

(2.7) ||^l||oo,||V^2||oo,||^;i|oo,||V^^||oo<C. 

The Wronskian of the solutions ifji and ip2 is independent of x: 

(2.8) W{^i,il)2}--=MxW2{x)-i^'i{x)il^2{x), for all .T e M+ \ X 

This is easy to check: it is constant on every interval (a;„,x„+i) and at the points 
{x„}rigN one has ''/'2}(2;n+) = 'i/'2}(a;n— ) from (|2.2p . The Wronskian 

is non-zero since the solutions are linear independent. From (|2.8p one has: 

IW'Wi, V'2}| < C(|^i(a;)| + for ah a; e M+ \ X, 

and therefore there exist constants C, and C* such that 

(2.9) G<C^<\iJi{x)\ + \i,[{x)\<C'\ for all a; e M+ \ X 

Now we apply the trick used in the proof of |S92[ Lemma 4]. We consider for 
an arbitrary n S N the interval [a;„,x„+i]. Since the function "01 is continuous on 
[xn, a;„4.i], the formula 

(2.10) 4>i{Xn+l) - MXn) = / 4>[{t)dt 

holds. Set 

Ax* 

p* 



Ax^ + 2 

Next we show that there exists a point x* G [x„,x„+i] such that |^/'i(x*)| > p^^C^:. 
Let us suppose that such a point does not exist, i.e. |V'i(2^)l < P*C'* for all x G 
[xn,Xn+i]- We get from (|2.9p that ji/'U^^)! > (1 —p*)C■^, for every x G [xn,Xn+i]- In 
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particular ip[ is sign-definite in [xn,Xn+i], so using (|2.10l) and Aa;„ > Ax^, we get 
a contradiction 

2p,C, > \'ilJi{xn+i)\ + \ipi{xn)\ > |V'i(a;„+i) - 'lpl{Xn)\ = 

\i)[{t)\dt > Ax„(l -p*)C^ > Aa;*(l -_p*)C, = 2p*C*. 



Thus the point a;* with required properties exists. 

Since IV'U^;)! < Ci for every x G [a;„,a;„+i] such that |a; — a;*| < one has 



P,C. 



We have shown that every interval [xn, Xn+i] contains a subinterval 



of length I := min (Ax*, 2f§^), where > Therefore 



\Mt)?dt> 



2 

2r'2 



On the other hand, 



IV'iWprfi < Ax*C'^n. 



Summing up, for every solution ij} the integral \t}j{t)\^dt has two-sided linear 



estimate. Thus no subordinate solution exists. 



□ 



2.4. Benzaid-Lutz theorems for discrete linear systems in C^. The results 
of |BL87j translate classical theorems due to N. Levinson [L48 and W. Harris 
and D. Lutz |HL75) on the asymptotic integration of ordinary differential linear 
systems to the case of discrete linear systems. The major advance of these methods 
that they allow to reduce under certain assumptions the asymptotic integration of 
general discrete linear systems to almost trivial asymptotic integration of diagonal 
discrete linear systems. For our applications it is sufficient to formulate Benzaid- 
Lutz theorems only for discrete linear systems in C^. The first lemma of this 
subsection is a direct consequence of |BL871 Theorem 3.3]. 

Lemma 2.3. Let /i± € C \ {0} he such that ^ and let {Vn}neJi € 

£^(N, C^^^). If the coefficient matrix of the discrete linear system 



It 



n+l 











is nan- degenerate for every n E then this system has a basis of solutions it^ 
with the following asymptotics: 



It- 



oil) 



]^(^+ + (14)ii) as n -)■ oo. 



k=l 



)22) as 



k=l 



where by (Vfe)ii and (^4)22 we denote the diagonal entries of the matrices Vk, and 
the factors {p,+ + (Vfc)ii) and + {Vk)22) should be replaced by 1 for those values 
of the index k for which they vanish. 



The following lemma is a simplification of |BL87[ Theorem 3.2]. 

Lemma 2.4. Let {t„}„eN G ^^(N) and {KjneN & C^^^) satisfy the following 
conditions. 

(i) Levinson condition: there exists M_ such that for every ni > 712 



n ii-t„i <M_ 
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and either there exists such that for every ni > ?i2 



n ii+i„i > 



or 



|1 + tfel oo as n ^ CO. 
(ii) The sum ^ Vn is conditionally convergent and 



If for every n G N 



det 



l + tr. 







1 - 1, 



Vn 



^0, 



then the discrete linear system 



l+t„ 







1-tr 



Vn 



~^n 



has a basis of solutions ~it^ with the following asymptotics: 



o{l] 



as n 



oo, 



k=l 
-t n 



0(1) 



as n -)■ OO, 



k=l 



where the factor (1 — tk) should be replaced by 1 for those values of the index k for 
which it vanishes. 

2.5. Characterization of spectra for operators with point interactions via 
Weyl function. In this section we view the operator H^ x,a as a self-adjoint 
extension of a certain symmetric operator and we parametrize this self-adjoint 
extension via suitable ordinary boundary triple |Ko75[ IBr76) . As a benefit of 
this parametrization we reduce the spectral analysis of the self-adjoint operator 
H>c,x,a to the spectral analysis of an operator-valued function in the auxiliary 
Hilbert space. The reader is referred to |DM95[ IBGP08) for an extensive dis- 
cussion of ordinary boundary triples and their applications. The Weyl function 
method was applied to Schrodinger operators with (5-interactions by many authors 
[KH891 IMMl IKMIOI iGOTOl lAKMlO) . 

Let X = {xn : n S N} be a discrete set satisfying the condition (|2.ip and let k 
be a parameter in the interval [0,7r). Define the minimal symmetric operator via 
its action and domain 



cos(>r)'(/'(0)— sin(x)'0 (0) 



L 4>'(x,^-)=ip'(x„ + ), nGN 

It can be easily verified that H^ x is a closed symmetric densely defined operator 
with the adjoint (the maximal operator) 



H* 



dom HI 



■ X 



'(0) 1 
GN. / 
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Definition 2.5. For the discrete set X and the parameter >c as above we define 
^o,>c,x and Ti^^^x as 



By [BGPOSl Proposition 3.17] the triple IV^^x {l'^{^),To,.<,x,Ti,x,x} is an 
ordinary boundary triple for H'^ ^ in the following sense: 



(i) ran 



(ii) for all f,g e dom H*^^ 



Remark 2.6. The assumption inf |a;„-|-i— a:„| > is essential. In a recent work IKMlOj 

nGN 

an ordinary boundary triple was constructed for a class of sets X not satisfying this 
condition. 

For further purposes we define the set cr^ x := '^{H^^ \ kerFo.^.x)- We asso- 
ciate the Weyl function with the triple li^^x 



(2.11) 



M^,x{\) ■.= Ti,^,x{Ta,^,x rker(7?: 



,x 



XI) 



which is defined for all A € C \ cr^ x- 

Recall that Sn{k) = sin(A:(a::„_(-i — Xn)) and c„(fc) = cos(/c(a;„+i — a;„)). Given a 



vector ^ e £ (N) we can find a solution rp e ker(i/^ ^ — XI) such that ■i/;(a 



The operator AI^_x{X) maps to the vector {ip'{xn + 0) — ip'{xn — 0)}neN & £'^{N). 
A straightforward calculation based on the formula (j2.4l) gives us the values of all 
entries of the matrix M^^xiX) except the first row. To find the values of the first 
row entries one needs to know 4'{x) on the interval [0,a;i], i.e., to solve the non- 
homogenious Sturm-Liouville problem —^"{x) = X4>{x), V'(O) cosx — ip^O) sin^r = 
0,ip{xi) = —£,1. The result of calculations is the following: 



M^,x(fc2) 



/foi(fc) ai{k) 

ai(fc) 62 (A:) 02 (/c) 

02 (A:) bsik) 

V ; ; ; 



ith 



biik) : 

bn{k) 

an{k) 



= k 



Co{k) cos >r — ksQ^k) sin >f 
so{k) cos >c + kca{k) sin >c 

I Cn{k) ^ C„_i(fc) \ 

ys„(fc) s„-i(fc) J ' 

k 



...\ 

. . . 1 

■■} 

Cljk) 

's,{k)' 



n > 2, 



s„(fc) ' 



Remark 2.7. Note that the set i7>c,jc consists of squares of all zeroes of denominators 

in the formulas above: 

(2.12) 

a^^x — {k^ G ]R : Sn{k) — for some n e N or so(fc) cos >c + fcco(fc) sin >c = 0}. 



£2(N) 

(2.13) 



Let a — {a„}„gN be a bounded sequence of real numbers. Define an operator in 

= {-anin)nen, e e ^'(N), 
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which is bounded and self-adjoint. The operator H^^x,a is a restriction of the 
operator H^^^ corresponding to in the fohowing sense 

H>c,x,a — H^,x \ ker(ri_^_x — BaTo^>t,x)- 

Proposition 2.8. |BGP08[ Theorem 3.3] Let the discrete set X , the sequence of 
strengths a and the self-adjoint operator H^t^x^a be defined as in Section \2.1\ Let 
ly^.x be the set defined in (j2.12[) and Ba be given by (|2.13p . Let Af^^.x be the Weyl 
function (12.111) associated with the ordinary boundary triple T\^ x- For A G M\cr>f^jf 
the following holds: 

X e acss{H^^x,a) if and only if e (Tcss{Ba - M^,x(A)). 

3. Spectral and asymptotic analysis 

3.1. Asymptotic analysis of a special class of discrete linear systems. In 

this section we study a special class of discrete linear systems that encapsulates 
system (|2.5p corresponding to X — {xn : n g N} and a = {anjneN as in Models I 
and II described in the introduction. 

Let the parameters /, a, b, 7 and uj satisfy conditions 

(3.1) leM.+ , a,beC, and a; e ( 0, 

For further purposes we define 
(3.2) 



(3.3) 
and 



fi± I ± y/P - 1, 
l3 |z|, 



2i sin w 



argz, 



(3.4) 



exp 



\ 1-7 



if 7 < 1, 

n^'", if 7 = 1. 

The following has technical nature and helps to simplify the analysis of cases (Model 
I and Model II). 



Lemma 3.1. Let the parameters I, a,b, 7 and uj be as in (j3.ip . Let fi±, /3,ip and 
fni ) be as in (j3.2l) . (13. 3p and (13.41) . respectively. Let the sequence of matrices 
{Rn}nefi G ^"'^(N, C^^^}. If the coefficient matrix of the discrete linear system 



(3.5) it 



n+l 





-1 



0\ e 
a b 



0\ e 
a h 



R„ 



n't \a bj n'l 

is non-degenerate for every n e N, then this system has a basis of solutions li^ 
with the following asymptotics. 

(i) Ifl€M.+ \ {coscj, 1}, then 

1 



0(1) 



as n ^ 00. 



(ii) If I = cosw, then 



It: 



cos(wrt + (^9/2) 
cos(w(7i + 1) + (^5/2) 

sin(wrt + (^5/2) 
sin(w(n + 1) + ip/2) 



-o(l) 



0(1) 



fniP) as n->- 00, 



fn W) asn-¥ 



Before providing the proof we make several remarks. 
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Remark 3.2. We do not consider the (double- root) case I — 1 . The analysis in this 
special case is technically involved and we refer the reader to |J06[ IJNSh07| . 

Remark 3.3. The case of negative values of the parameter I can be treated by the 
same lemma if one makes the substitution 



This leads to the same answer in the case I e R_ \{ — 1, — cosw} and to the following 
answer in the case / = — cosu: 



and 

where 
(3.6) 



cos(wn + 
— cos(ti;(n + 1) + 



sin(a;n + (^1/2) 
— sin(w(n + 1) + <(5i/2) 



-o(l) 



0(1) 



/+(/3i) asn^oo, 



be 



-2iuj 



zi 



2i sin oj 

The (double-root) case Z = — 1 is out of our considerations. 



argzi. 



Proof of Lemma \3.1\ Since I e \ {1}, the constant term in the coefficient matrix 
in (j3.5l) can be diagonalized as follows 



1 1 
In view of the identity 



1 

-1 21 



1 



1 



^l+ 
^_ 



1 1 



the substitution 
(3.7) 



lin = 



1 1 

transfornis the system (j3.5p on lt„ into the system on if „ given below 

'-(a + /i+6) -(a + Ai-fe)\ e"'"" 



It 



n+l 



(3.8) 



M+ 




1 



-(a + /i+5) -(a + M-&)A 



(i) The case I E IR+ \ {cosw, 1} splits into two subcases: / > 1 and / G [0, 1) \ 
{coso;}. 

The condition I > 1 implies < ^_ < 1 < fi+, thus Lemma [^751 is applicable to 
the system (|3.8|) and it gives us a basis. Reverting the substitution p.7p we get the 
statement. 

The condition I e [0,1) \ {cosw} implies j/i+j = = 1. We reduce this 

situation to Lemma 12.41 by the following substitution: 



(3.9) 



^ll 
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The system on lit„ has the foUowing form 
(3.10) 



1 



A*- - AH VM+(a + Ai+^)A'4 



Now Lemma HH] is apphcable and reverting substitutions p.7p and p.9p we get the 
statement. 

(ii) We start with the system p.lOp . In the case / = cos uj Lemma 12.41 could 
not be apphed immediately. In order to bring the system into the form in which 
Lemma 12.41 is applicable we make further transformations. Since I — cosw, we, 
clearly, have fj,± = e^*" and thus, grouping summands in a right way, we can 
rewrite our system in the form 



U^n+l 



Z j „7 + 



with 
and with 



be 



-2iuj 



2i sin oj 



1 



(a*- - 



2iijn 



-^-{a + ii+b) -^_(a + 

^+(a + /i_6) 



R„ 



The sequence {Vn\n&i satisfies the condition (ii) in Lemma 12.41 It can checked by 
a straightforward calculation that |z| = /3 and arg(z) — with /3 and Lp defined in 



Since 

the substitution 
(3.11) 



(ie'^P 
Pe^^'fi 



/3 
-/3 



leads to the system on ~^ n of the form 



(3.12) 



n+l 



P 

-/3 / ,^7 



with {KilngN satisfying the condition (ii) in Lemma 
For sufficiently large no 



(3.13) 



n ^^^^ ' -^"^^^ as n ^ cx), 

k=no ^ ^ 
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with fn{-) defined in (13.41) . Now we apply Lemma [2.41 to the system p. 121) and 
using p.l3p we get a basis of solutions of that system of the form 



(3.14) 



0(1) 



0(1) 



/+(/3) asn^oo, 



/„ (/?) as 



n — )■ oo. 



The asymptotics above imply the statement immediately after reverting substitu- 
tions (|3.7p . p.9p . (13. lip and coming back to the sequence itn- D 

3.2. Decomposition of the transfer matrices. Recall that Tn(k) denotes the 
transfer matrix (j2.6l) . In this section we decompose this transfer matrix subject to 
Model I and Model II in the form 



(3.15) 








a b 




a b 



where i?„ £ £^(N;C^^^). This decomposition allows to apply later Lemma [01 

Model I: Amplitude perturbation. The sequence of strengths a = {a„}„gN and 
the discrete set X — {xn'- n G N} are defined as follows 



csin(2w7i) 
ao H \- q-a, 



Xn :— nd and x„+i — Xn = d. 



Substitution of these expression into (j2.6p and tedious calculation show that T„{k) 
for k ^ nZ/d has the representation of the type (I3.15P with the following values of 
the parameters: 

ao sm{kd) 



(3.16) 



cos{kd) 



0, 



b = 



2k 
csin(fcc?) 
2ik ' 



Ls{k), 



Model II: Positional perturbation. The sequence of strengths a = {S„}„gN and 
the discrete set X — {xn : n E N} are defined as follows 

csin(2a;n) 



nd - 



nf 



In this case substitution into (j2.6p and calculation give for k ^ nZ/d the represen- 
tation of the type (j3.15l) with the following values of the parameters: 



I = cos{kd) 



aQ sin(/sd) 
2k 



= Lsik), 



(3.17) 



: — 2ikc cot{kd) sin^ (w) , 

-c cos{kd) sin(a;) 4k cos(cj) cot{2kd) — 2k cot{kd)e' 



Remark 3.4. We do not provide here long calculations for the decompositions of 
the transfer matrices. In these calculations we in fact expand the transfer matrix 
in powers of -^^ as n — >■ cxd up to a summable term. 

3.3. Asymptotics of generalized eigenvectors for the spectral bands. In 

this section we obtain using Lemma |3T] asymptotics of solutions ^ of the difference 
equation (12.31) in the cases of Model I and Model II. Using the statements from 
Section 12.31 we come from the asymptotics of these solutions to the conclusions 
about the spectral bands and critical points inside them. 
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In what follows we use the inverse Joukowsky mapping $ : C i— t- C defined as 
(3.18) :=w+ Vu'^-l, weC, 

where the branch of the square root is chosen so that > 1 for w e C+. In 

the table below we list certain functions playing a role in the main theorem of this 
section. 



Table 1 . Notation for Theorem 13.51 



Perturb.: 

Funct.: 


Amplitude 


Positional 


m ■■= 


c sh'i{kd) 
4:k sin LJ 


1 cao 1 
1 2 1 


d±{k) := 


iarg(T^) 


i arg(icao) 



Theorem 3.5. Let the parameters d, a^, c, uj, 7 and {qn}nen (jl.2l) . Let 

the function Lg be as in (jl.7p . Assume either that X — X, a = a as in Model I or 
that X = X , a — a as in Model IL Then for A G R the finite difference equation 
(12. 3p with k — a/a has a basis of solutions £,^{k) — {C^(fc)}„gN with the following 
asymptotics: 

(i) If\Ls{k)\ e M+ \ {cosw,!}, then 

e,t(fc) = ($(£A-(fc)))^"(l + o(l)) as n ^ 00, 
where $ is defined in (|3.18p . 

(ii) If Ls{k) =coscli, then 

^+(fc) = (cos(^Jr^ + ^?+(fc)) +o(l))/+(/3(fc)) as n 00, 

C{k)^{sm{un + d+{k)) + o{l))f-{m) as n ^ ^. 

If Ls{k) — — cosw, then 

^+{k) = (-!)«( cos(c.n + + o(l))/+(/3(fc)) as n^^, 

C(fc) = (-l)"(sin(c.n + + o(l))/-(/3(fc)) as n^^. 

Here f^i') defined in (j3.4l) and the functions /?(•) and t?±(-) are given 
in the 1st column of the Table 1 in the case of Model I and in the 2nd 
column in the case of Model II, respectively. 

Proof. Let us fix fc = \/X such that 7^ 1. For such k the limit property 

lim inf „_i.+oo > holds, and the transfer matrix Tn{k) of the discrete linear 

system \2.b\ decomposes into the form in (j3.5|) with the parameters a, b and I given 
in (|3.16p for the Model I and given in (j3.17p for the Model II, respectively. 

In the case S K+ \ {coso;,!} Lemma IXTl fi) could be applied and it 

automatically implies the item (i). Assume now Ls{k) = zbcoscj. In the case of 
Model I we plug into formulas p.3p and p. 61) the values a and b from p.l6p and 
we easily get the following 
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In the case of Model II we plug into formulas p.3p and p. 61) the values a and b 
from (|3.17p and usmg that Ls{k) — ± coscj after a long and tedious calculation we 
get 

z{k) = zi(fc) = . 

Now Lemma I^TTI fii') could be applied. It gives us the asymptotics for the sequences 
of C^-vectors 



ltt{k) = 



Extracting the second components from these asymptotics we obtain the statement. 

□ 

Given the asymptotics, we can now come to the conclusions about the structure 
of the purely absolutely continuous spectrum using the subordinacy theory, see 
Section [^31 Let the parameters d, ao,c, co and 7 be as in (|1.2p . Recall that the self- 
adjoint operator corresponds to the Kronig-Penney without perturbation, where 
the distance between interaction is the constant d and the strength of interactions 
is the constant ao- The spectral properties of this operator are discussed in the 
introduction. The absolutely continuous spectrum of the operator is the set 

'J.ciH^) := {agM: L5 (Va) G [-1,1]}. 

Denote the set of all critical points by 

Scr {a e M: Ls (Vx^ = icoswj C Int(CTac = |a e M: Ls (Vx^ G (-1,1)}- 

Corollary 3.6. Let the parameters d, ao,c, uj, 7 and {<?„} be as in (jl.2l) and let 
the sequences X,a and X,a be as in (|1.3p . (|1.4p and in p.5p . (|1.6p . respectively. 
Then the following statements hold. 

(i) The spectrum of the operators ^ ~ and s purely absolutely con- 
tinuous on the set 

Int(aac(i?>.)) \ 

(ii) Let A = fc^ G Scr- If either 7 < 1 or 



7 = 1 and 



csm(kd) 



2k sine 



> 1, 



then there exists a unique value k such that A is an embedded eigenvalue of 



H 

(iii) Let A = fc^ G Scr- if either 7 < 1 or 

7 = 1 and |cao| > 1, 
then there exists a unique value >c such that A is an embedded eigenvalue of 

Proof. For fc^ G Ini {a ac{H>cS) one has that liminf„^oo |'Sn(fc)| > 0. From the 
asymptotics of (^) ^-nd the formula (12. 4p it follows that all the solutions of the 
spectral equations for both operators ^ ~ and ^ ~ are bounded. The asser- 
tion of the item (i) follows now from Proposition 12.11 Proposition 12.21 and bound- 
edness of all solutions. 

Conditions of the items (ii) and (iii) guarantee the existence of the square- 
integrable solution of the spectral equations at the point A for the operators ^ ~ 
and ^ g, respectively. This fact means that there is a unique value of the bound- 
ary parameter x such that this solution satisfies the boundary condition, and thus 
in this case A is an eigenvalue. 
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□ 

3.4. Spectral gaps. In this section we show that the essential spectra of perturbed 
operators coincide with the essential spectrum of the non-perturbed Kronig-Penney 
model. Let the sets ^ and ^ be defined by X and X, respectively, by the 
formula (|2.12p . Recall that the Weyl functions ^ and ^ defined in 
(|2.1ip and operators Bs and Bs are defined in (|2.13p . 

Theorem 3.7. Let the parameters d, aQ,c, lo, 7 and {q,i\ be as in (|1.2p . Let the 
sets X,X and the sequences a, a be as in (jl.3p . (|1.5p and (|1.4p . (|1.6p . respectively. 
Let the self-adjoint operator correspond to the non-perturbed Kronig-Penney 
model. Let the operators ^ ~ and ^ ~ correspond to X,a and X,a as in 
Section \2.1l respectively. Then the spectrum of both operators ^ ~ and 2 s 
on the set M \ acss{H^) is pure point. 

Proof, (i) Define Bag := —aol, a diagonal operator in ^^(N). One has: an cto 
as n — 00. Thus 

Ba ^ ®oo- 

Consider A e M\ {a^ss{H^)Ua^j^) . By Proposition EJl i a^ssiBao - M^j^{X)) . 
Therefore by compact perturbation argument ^ acss{Bs — ^^^{X)) and again 
by Proposition 12.81 A ^ acss{H^ ^ ^. Thus the spectrum of the operator ^ ~ 
on the set K \ ((Tcss(-ffx) U j^) is discrete. Since the set cr^ n (R \ crcss{H^)) is 
at most countable, it cannot support absolutely continuous or singular continuous 
spectral measure. 

(ii) One has Bs — Bag . Consider A e ]R\ (crcss(-ff>t) U U cr^ j-= ) . Let 5„(fc) be 
diagonal and a,n{k) be anti-diagonal entries of the Jacobi matrix x C-^) (where 
A = k^). Clearly a„(fc) — > ^-^^^^^ and bn{k) — > 2kcot{kd) as n — > 00. It follows 
that 

Again, Proposition 12.81 implies that A ^ Ccss {H^ ^ s) • Hence the spectrum of 
^x,x,s on M \ (acssiH^) U a^j^ U a^j^) is discrete. The set (cr^j^ U a^j^) n (M \ 
aess(H^)) is at most countable, so the spectrum of ^ g on the set M \ acss{H^) 
is pure point. □ 
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